Nonequilibrium carriers in an intrinsic graphene under interband photoexcitation 
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We study nonequilibrium carriers (electrons and holes) in an intrinsic graphene at low temper- 
atures under far- and mid-infrared (IR) radiation in a wide range of its intensities. The energy 
distributions of carriers are calculated using a quasiclassic kinetic equation which accounts for the 
energy relaxation due to acoustic phonons and the radiative generation-recombination processes 
associated with thermal radiation and the carrier photoexcitation by incident radiation. It is found 
that the nonequilibrium distributions are determined by an interplay between weak energy relax- 
ation on acoustic phonons and generation-recombination processes as well as by the effect of pumping 
saturation. Due to the effect of saturation, the carrier distribution functions can exhibit plateaus 
around the pumping region at elevated intensities. As shown, at sufficiently strong mid-IR pumping, 
the population inversion can occur below the pumping energy. The graphene dc conductivity as a 
function of the pumping intensity exhibits a pronounced nonlinearity with a sub-linear region at 
fairly low intensities and a saturation at a strong pumping. However, an increase in the pump- 
ing intensity in very wide range leads only to a modest increase in the carrier concentration and, 
particularly, the dc conductivity. The graphene conductivity at mid-IR irradiation exhibit strong 
sensitivity to mechanisms of carrier momentum relaxation. 
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PACS numbers: 73.50.Pz, 73.63.-b, 81.05.Uw 



I. INTRODUCTION 



The features of the dynamics of carriers (electrons and holes) in graphene and the mechanisms of their relaxation 
[l| result in the exceptional properties of graphene and a new device prospects (see 0] for review). The studies of 
optical phenomena, including the Raman scattering (see0] and references therein), ultrafast spectroscopy, 0,0] and 
magnetooptics, Q can be used to reveal both the energy spectrum parameters and the mechanisms of carrier scattering. 
The gapless energy spectrum of graphene, with the characteristic velocity vyy — 10 8 cm/s which corresponds to 
the neutrinolike bandstructure (Weyl- Wallace model) Q, provides its nontrivial optical properties due to effective 
interband transitions in the far- and mid-infrared (IR) spectral regions. The linear response of epitaxial graphene Q 
and graphite 0] was measured and the pertinent calculations were performed in Ref. 10 1. Recently, 11 1 a fairly low 



^ ■ threashold of the nonlinear response under far- or mid-IR excitation have been found and the photoconductivity of an 
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intrinsic graphene was calculated for a low-pumping region. But an essentially nonlinear regime of response was not 
calculated and no experimental data concerning nonlinear properties of graphene in this spectral region are available. 
- 1 Thus, an investigation of the nonlinear response of graphene is timely now. 

In this paper, we consider the nonequilibrium energy distributions of the carriers under interband photoexcitation by 
far- or mid-IR pumping in a wide range of the intensities. Using the obtained distributions, we analyze the nonlinear 
dependencies of the carrier concentration and the dc conductivity on the pumping intensity and frequency at different 
temperatures. The dynamic conductivity of graphene photoexcited by far- or mid-IR radiation is calculated as well. 

As shown below, the energy distribution of carriers are determined by an interplay between quasielastic energy 



'Electronic address: ftvasko@yahoo.com 



2 



relaxation, which is ineffective at low energies, and generation-recombination processes. There is a marked increase in 
the carrier population in the range of energies where the phonon and radiative mechanisms of relaxation are virtually 
compensated. The interband absorption saturation also affects the energy distribution of carriers. The graphene dc 
conductivity as a function of the pumping intensity exhibits a pronounced nonlinearity with a sub-linear region at 
fairly low intensities and a saturation at a strong pumping. It is also found that in a certain energy range of the 
pumping power the real part of the interband contribution to the dynamic conductivity becomes negative, i.e. the 
negative absorption takes place. 

The paper is organized in the following way. The model under consideration is described in Sec. II. In Sec. Ill, wc 
present the nonequilibrium carrier distribution as function of pumping intensity for a few pumping frequencies. The 
results of calculations of the photoconductivity and the dynamic conductivity, as a response on a weak dc electric field 
or a probe high-frequency field, are discussed in Sec. IV. The brief discussion of the assumptions used in calculations 
and concluding remarks are given in the last section. 



II. MODEL 



In order to obtain the energy distributions of the carriers in intrinsic graphene we use the quasiclassic kinetic 
equation derived and analytically analyzed previously for the case of low intensity pumping. The kinetic equation 
under consideration accounts for the energy relaxation due to scattering on acoustic phonons, the radiative generation- 
recombination processes associated with thermal radiation, and the far- or mid-IR pumping. We disregard the inter- 
carrier scattering since at low temperatures the carrier concentration can be small even at relatively strong interband 
pumping. Since the scattering mechanisms in c- and w-bands are symmetric, the electron and hole distributions in 
the intrinsic material are identical and we consider below the carrier distribution function f p . Taking into account 
the abovementioned mechanisms, the kinetic equation under consideration, which governs the distribution function 
f p , is presented in the following form: 



JlaW} + Mfp} + G{f p } = 0. 



(1) 



Here the collision integrals Jla{Ip} and Jr{/ p } are associated with the relaxation of carriers caused by the acoustic 
phonons and the equilibrium thermal radiation, respectively, the term G{f p } describes the interband carrier excitation. 
Since the interband transition due to the acoustic phonon scattering are forbidden (the sound velocity is weak in 
comparison to vw), the concentration balance equation takes form 



d PP [J R {fp} + G{f p }} = 



(2) 



and it can be considered as the normalization condition for f p . Another condition for f p streams from the requirement 
so that any term in Eq. (1) is equal to zero at p — > oo, i.e. the zero acoustic flow at high energies takes 
place, Jla{/ p }| p ->oo=0. 

Considering the quasielastic scattering of carriers on acoustic phonons, one can use the Fokker-Planck form of 
Jla{Ip\ [12j and Eq. (1) can be presented as 11 1 
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Here pr = T/vw is the characteristic thermal momentum (T is the temperature in the energy units), iV 2p / PT = 
[exp(2p/px) — 1] is the Planck distribution function, and the form-factor A(e) describes the broadening of interband 
transitions. So that pn — hQ/2vw is the momentum of just photogenerated carriers corresponding to the pumping 



3 



frequency fl and the broadening is described by 5pn which is proportional to pa if the scattering rates are proportional 



to the density of states. It was assumed below that A(e) = tt ' exp(— e ) and Spa/pn —0.1 [13J. The dimensionless 



parameter 7 and G are the relative strength of the energy relaxation with respect to the generation-recombination 
efficiency and the characteristic value of the pumping intensity respectively. It is important below, that a following 
dependencies on T, Q and the pumping intensity S, take place: 



7 oc T 



G(xS/(Q 3 T), 



(4) 



so that the acoustic contribution increases with T while the excitation efficiency decreases with and T. The explicit 
expressions of g and G can be found in Ref. Using the typical parameters of graphene [l4[, at T=77 K one 

obtains 7 ~0.32. At the same temperature at S—YW/cm 1 and 7j,O=120 meV one can get G ~2.8. 

Numerical procedure for the problem described can be simplified if one takes into account that the equilibrium 
distribution remains valid at p — > 0. The condition f p ^ Q = l/2 (it is the Fermi function at zero energy) can be explicitly 
derived considering that G{f p } vanishes at p — > (if pn ^> Spn) and equations Jla{/p} — and </r{/ p } — are 
satisfied by the equilibrium distribution. Taking this into account, one can write the boundary condition at p — > 00, 
as the requirement that the factor {. . .} under derivative in the first term in the left-hand side of Eq. (3) turns to 
zero. As a consequence, for numerical solution of Eq. (2), one can use the following requirements as the boundary 
conditions for the distribution function: 
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(5) 



while the balance condition (2) should be used for checking of a numerical results obtained. The numerical solution 
of Eq. (3) is performed below using a finite difference method and the pertinent iteration procedure. 
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FIG. 1: (Color online) Distribution functions f p vs energy vwP (a) for the pumping intensities 5 = 10 



10" 



10 3 W/cm 2 



(from left to right) at T = 4.2 K and (b) for 5=10~ 2 , 10~\ . . ., 10 3 W/cm 2 at 77 K. Dotted curves correspond to the equilibrium 
distributions. Vertical arrow corresponds to the energy equal to T1Q./2 at hQ.=12 meV. 



III. NONEQUILIBRIUM DISTRIBUTION 

We present here the results of numerical solution of Eq. (3) with conditions (5) and discuss the obtained distribution 
functions at different excitation conditions (frequency and intensity of pumping) and temperature for the typical 



parameters of graphene 1J]. The variation of the sheet carrier concentration with varying excitation conditions is 
also considered. 

The obtained distribution functions under the far-IR pumping with the photon energy Ml =12 meV are shown in 
Fig. 1 with the one-order step in pumping intensities from 1CU 5 W/cm 2 to 10 3 W/cm 2 . First of all, one can see 
a visible modification of the carrier distributions at low intensities (up to S ~ 10~ 3 W/cm 2 at T—4.2 K and for 
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FIG. 2: (Color online) The same as in Fig. 1 under mid-IR pumping with TiQ. =60 meV. 
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FIG. 3: (Color online) The same as in Fig. 1 under the CO2 laser pumping, TiQ. =120 meV, with intensities (a) 5=10 
10" 3 , . . ., 10 5 W/cm 2 (from left to right) and (b) for S=10 -2 , HT 1 , . . ., 10 5 W/cm 2 . 



S < 10 _1 W/cm 2 at T—77 K). Next, the peak of the distribution function shifts towards Ml/2. Such a behavior is in 
agreement with the previous analytical consideration Further, at sufficiently strong pumping, f p ^ pn tends also 
to the values close to 1/2, so that plateau-like energy distibutions are formed. The range of carrier energies, where 
f p ~ 1/2, widens with increasing pumping intensity starting S ~ 10~ 2 W/cm 2 at T—A.2 K and S ~ 0.1 W/cm 2 
at T=77 K. Since f p remains equilibrium at p — > (all terms of Eq.(l) are equal zero separately), one can see a 
deepening of f p in the region vwP ~ T. 

Under mid-IR excitation, hil =60 meV, a similar character of the low-pumping regime of response takes place, see 
Fig. 2. With increased pumping intensity, a peak of distribution appears, moreover, the peak distribution function 
can markedly exceed 1/2, i.e. the population inversion occurs if S > 10~ 2 W/cm 2 at T=4.2 K and S > 10 _1 W/cm 2 
at T=77 K. Further, in the pumping region when S > 1 W/cm 2 at T=4.2 K (or S > 10 W/cm 2 at T=77 K) f p ^ pn 
tends to 1/2, so that a plateau of f p is formed around the energy Ml/2. As S increases and the plateau region widens, 
the peak amplitude below hQ/2 somewhat decreases but a maximal value of f p exceeds 1/2. Since the equilibrium 
distribution is predetermined by Eqs. (3) and (5) at p — > 0, a non-monotonic distribution with a deepening at 
vwP < T, a peak under the energy Ml/2, and a plateau around hQ/2 is realized. 

Similar character of distribution takes place as Titt increases: the same peculiarities are shifted to higher intensities 
according to Eq. (5). In Figs. 3a and 3b we plot f p under the CO2 laser pumping, Ml =120 meV. One can see 
that the population inversion regime begins for S > 10~ 2 W/cm 2 at T=4.2 K (or for S > 10 _1 W/cm 2 at T=77 K) 
and the saturation region around Tiil/2 takes place if S > 10 3 W/cm 2 at T=4.2 K and 77 K. These peculiarities are 



5 



retained for higher intensities, up to S ~ 10 5 W/cm 2 . 
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FIG. 4: (Color online) Carrier concentration n2D vs pumping intensity S for the excitation energies hfl =12, 60, and 120 meV 
(solid, dashed, and dotted curves, respectively) at different temperatures. 

The obtained distribution functions allows to calculate the sheet carrier concentration at different pumping condi- 
tions and temperatures according to the standard formula: 

2 f°° 

n 2 D = —2 / dppf p . (6) 
irn Jo 

Figure 4 demonstrates the dependences of the sheet concentration in graphene as a function of the pumping intensity 
corresponding to the energy distributions shown in Figs. 1-3. As follows from the numerical calculation of the integral 
(6) at T=4.2 K, n-iB increases with S in an intermediate pumping region and tends to saturation under higher pumping 
intensities. The power dependence n2D oc S r with r ~0.57 takes place at T = 4.2K between S ~ 10~ 6 W/cm 2 and 0.5 
W/cm 2 for the mid-lR pumping or between S ~ 10~ 7 W/cm 2 and 10 -5 W/cm 2 for the far-lR pumping (see Figs. 4a 
and 4b, respectively). Under mid-lR pumping at T=77 K, the concentration increases according to the same power 
law and saturates at the same pumping intensities. However, the nonlinear regime begins starting S ~ 5 mW/cm 2 . 
This is because the equilibrium concentration is proportional to T 2 . Note that in the case of pumping by CO2 laser, 
the saturation occurs at n2D > fO 11 cm~ 2 . In this situation, the inter-carrier scattering might be important. 



IV. DC AND DYNAMIC CONDUCTIVITIES 



Here we turn to the consideration of the response of the nonequilibrium carriers with the distribution functions 
obtained above to a weak dc electric field or a probe radiation. Taking into account that the momentum relaxation of 
carriers is caused by elastic scattering mechanisms, one can use the following formula for the dc conductivity a 0: 



a = a 



2/ P 



dpf p 



*'(pJc/ft) 



(7) 



Here l c is the correlation length of static disorder scattering, 9(z) = e~ z2 Ii(z 2 )/ z 2 , where Ii(z 2 ) the first order 
Bessel function of imaginary argument, and Co is the conductivity in the case of short-range disorder scattering, when 
lc — 15j. According to Eq.(5) f p =o = 1/2, for the short-range scattering case, when pl c /fi <C 1 (p is the average 



momentum), the dc conductivity a is determined by the low energy carriers so that a appears to be independent 
of optical pumping intensity despite a significant concentration of the photogenerated carriers, a ~ Co- For the 
definiteness, it was assumed that l c — 10, 20, and 30 nm. 
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FIG. 5: (Color online) Normalized dc conductivity a/a vs pumping intensity 5*: (a) for low pumping region and Z c =20 nm at 
different temperatures and hQ = 12, 60, and 120 meV (solid, dashed, and dotted curves, respectively), (b) for mid-IR pumping 
at temperatures T = 4.2 K and 77 K (solid and dotted curves), and (c) for different l c and TiQ. =60 meV at T — 4.2 K and 
77 K (solid and dotted curves). 



Figure 5 shows the dependences of the dc conductivity a normalized by its characteristic value Co on the pumping 
power S calculated using Eq. (7) with the distribution functions obtained in Sec. Ill for T — 4.2 K and 77 K 
and different pumping energies Ml. The dependences are markedly nonlinear beginning from rather low pumping 
intensities (see the pertinent curves in Fig. 5a). The nonlinearity of the dc conductivity at low S (<1 mW/cm 2 ) 
is particularly strong in the case T = 4.2 K and far-IR pumping, Ml — 12 meV, but the variation of conductivity 
does not exceed 5% and becomes saturated at S >5 mW/cm 2 . In the range of intermediate intensities of mid-IR 
pumping, see Fig. 5b, there is a power increase in the dc conductivity, which is followed by the saturation region. It 
is instructing that an increase in the pumping intensity by several orders of magnitude leads to a modest increase in 
a/a . 

In Fig. 5c, the normalized dc conductivity as a function of the pumping intensity calculated for different values of 
the correlation length of the disorder scattering, l c . One can see that the conductivity becomes more sensitive to the 
pumping when parameter l c increases. If l c tends to zero, a quenching of photoconductivity takes place. Since the dc 
conductivity a depends on pl c /fi, it, as a function of the pumping intensity, saturatates at higher intensities that the 
carrier concentration n2D (compare Figs. 5 and 4). 

As seen from Figs. 2 and 3, the distribution function f p can be larger than 1/2 in a certain energy range if the 
pumping power is sufficiently high. This corresponds to the population inversion of carriers and might be used for far 
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FIG. 6: (Color online) (a) Spectral dependences of the real part of the dynamic interband conductivity Rea^ at T = 4.2 K for 
hfl=60 meV and pumping intensities: S = 10 -3 , 10 _1 , 10 1 , and 10 3 W/cm 2 (solid, dashed, dotted and dash-dotted curves, 
respectively), (b) Regions of parameters (shaded by vertical and horizontal lines for T—4.2 K and 77 K, respectively) for which 
the dynamic interband conductivity is negative if 7if2=60 meV. 



infrared lasing. The real part of the interband contribution to the dynamic conductivity which is given by Q 



= 2„,2 



P_ 

h 2 LU Jo \ ^ 



Rea. = I dpp(l - 2f p ) A ( Zj-ZZ ) , (8 ) 



can be negative. Here p^ = hu/2vw and Sp u oc p u are introduced in a similar way to pn and Spn in Eq. (3). 13J. At 
Sp^ — >0 one obtains Kea u < 0, if f Pul > 1/2. Figure 6a shows the dependences of the real part of the interband dynamic 
conductivity Recr w (normalized by factor e 2 /irh) on the probe frequency u at different pumping intensities S. In 
contrast to the pumping scheme when the pumping energy TiVL exceeds the energy of optical photons, so that a cascade 
of optical phonons is emitted and the photogenerated electrons (holes) are accumulated near the conduction band 
bottom (valence band top), the negative dynamic interband conductivity in the case under consideration corresponds 
to the interband transitions with relatively high energies hm (compare with 0). Figure 6 shows the regions (shaded) 
on the pumping intensity - probe frequency plane corresponding to the negative dynamic interband conductivity at 
different temperatures. The dependence of the real part of net dynamic conductivity on the probe frequency u>, can 
be markedly modified by the contribution of the intraband transitions (which correspond to the Drude conductivity). 
However, the latter contribution can be effective only in the range of rather small frequencies w.[l6| 



V. DISCUSSION AND CONCLUSIONS 



The main restriction the model used above is the neglect of the Coulomb interaction of carriers, so that the obtained 
results are valid at their sufficiently low concentrations. As can be found from Fig. 4, in the case of an intrinsic 
graphene, this assumption limits the validity of our model by the temperatures T < 100 — 200 K, whereas the limitation 
imposed on the pumping intensity is fairly liberal. The latter is owing to the effect of saturation of the interband 
absorption associated with the Pauli principle. In our consideration we also disregarded more complex processes, 
namely, the carrier interactions with the substrate vibrations and the effect of disorder on the generation-recombination 
processes. The mechanism included in the model describe the general features of the relaxation processes considered, 
so that their refinement should not lead to a qualitative change of the pattern of the phenomena studied above. As for 
the photogeneration of carriers by incident far- or mid-infrared radiation, only the single-photon interband processes 
were taken into account. Hence, possible nonlinear frequency multiplication and renormalization of the energy spectra 
were disregarded. Our estimates show that the latter correspond to the limitations of the pumping intensity by the 
values S < 0.5 MW/cm 2 (at HQ = 12 meV) and 5 MW/cm 2 (at HQ = 120 meV). The other assumptions of our 
model (the isotropic energy spectra of carriers, the valley degeneration, and exclusion of the interaction with optical 
phonons in the range of energy under consideration) appears to be rather natural. 
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In summary, we calculated the energy distributions of carriers in an intrinsic graphene at low temperatures under 
far- and mid-IR (frtt = 12 — 120 meV) pumping in a wide range of its intensities (from zero to 10 2 kW/cm 2 ). It 
was shown that these distributions are determined by interplay between weak energy relaxation on acoustic phonons 
and radiative generation-recombination processes as well as by the effect of pumping saturation due to the Pauli 
principle. The obtained energy distributions at elevated pumping intensities can exhibit the plateau regions and 
the regions corresponding to the population inversion. The graphene dc conductivity as a function of the pumping 
intensity demonstrates a pronounced nonlincarity at fairly weak pumping (particularly at helium temperature), a 
steep increase in a certain range of intermediate intensities, and a saturation at sufficiently strong pumping. Due to 
this, the graphene dc conductivity varies through 1-2 orders of magnitude when the pumping intensity is varied many 
orders. We showed that the alteration of the graphene dc conductivity with far- and mid-IR irradiation, i.e., the effect 
of the photoconductivity is sensitive to the correlation length of disorder scattering. 

To conclude, the obtained results, can be useful for the extraction of the parameters determined the relaxation 



mechanisms in graphene and for the estimation of a potential of novel optoelectronic devices [16l.ll7l.ll8l]. in particular, 
in terahertz and far infrared lasers [16L Il7f ■ 
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